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1 .  Introduction 


Goldfeld  and  Quandt  (1972)  consid¬ 
ered  a  switching  regression  model  with 
two  regimes: 

*i  “  ai  +  bi*i  +  uii 

yi  =  a2  +  b2xi  +  u2i 
2 

uli  ■v  IN  (o,  )  and  u2i  "v  IN  (o,o2*) 

They  assume  that  there  is  an  identifiable 
variable  Z^  such  that  if  Z^  <  ZQ  then  the 

observations  are  generated  by  (1)  and  if 
Z^  >  ZQ  then  the  observations  are  genera¬ 
ted  by  (2) .  They  then  define  the  step 
function 

D  (Z^  =1  if  Z.  <  Z. 

1  0  (3) 

=0  if  z.>z0 

Later,  instead  of  considering  the  step 
function,  they  consider  D(Z^)  to  be  a 
cumulative  normal  given  by 

D(Zi)  -  1  if  ZA  <  ZQ  +  ei 

*  0  if  zi  Z.  z0  + 

where  IN  (0,1). 

The  model  considered  by  Goldfeld  and 
Quandt  is  an  exogenous  switching  model 
because  is  assumed  to  be  independent 

of  Uj^  and  u2^.  Goldfeld  and  Quandt 

(1973)  give  an  extension  of  the  switch¬ 
ing  regression  model  to  simultaneous 
equations  systems  but  the  extension  is 
still  in  the  framework  of  exogenous 
switching. 

Maddala  and  Nelsocf  (1975)  discuss 
extensions  of  the  Goldfeld-Quandt  switch- 
regression  model  to  the  case  of  endogen¬ 
ous  switching.1  Many  of  the  practical 
problems  one  encounters  -  disequilibrium 
models,  selectivity  models  etc.,  are  all 
switching  regression  models  with  endoge¬ 
nous  switching.  The  present  paper  ex¬ 
tends  the  0-method  suggested  by  Goldfeld 
and  Quandt  to  switching  simultaneous 
equations  systems  with  endogenous  switch¬ 
ing. 
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There  have  been  previous  discus¬ 
sions  of  structural  change  by  Heckman 
(1976b)  in  models  with  endogenous  switch¬ 
ing.  However,  the  type  of  model  we  dis¬ 
cuss  has  not  been  considered  by  him. 

What  he  considers  is  a  shift  in  the  con¬ 
stant  term  (a  structural  shift)  and  not 
shifts  in  all  the  coefficients  as  we  do. 

What  the  D-method  does  is  to  com¬ 
bine  the  equations  in  the  two  regimes  in¬ 
to  a  single  equation.  This  will  enable 
us  to  apply  tests  for  equality  of  speci¬ 
fic  coefficients  in  the  two  regimes.  The 
alternative  to  this  procedure  would  be  to 
use  a  likelihood  ratio  test  every  time 
the  significance  of  a  particular  coeffi¬ 
cient  needs  to  be  tested.  Thus,  the  D- 
method  has  the  advantage  of  computational 
simplicity.  In  this  paper  we  suggest 
instrumental  variable  methods  to  estimate 
the  parameters  of  the  combined  equation. 
Following  the  methods  used  by  Amemiya 
(1978) ,  we  derive  the  asymptotic  covari¬ 
ance  matrix  (in  the  appendix)  and  present 
an  empirical  example  illustrating  the  use 
of  the  suggested  method. 

2.  The  Switching  Simultaneous  Equation 
System 

The  general  model  we  consider  in 
this  paper  is  the  following: 


I.  =  x!6  - 

E* 

1 

(4) 

Regime  1 : 

BiYii 

+  riZli 

=  eu 

if  Ii  >  0 

(5) 

Regime  2: 

B2Y2i 

+  r2Z2i 

=  e2i 

if  <  0 

(6) 

The  residuals  e^,  e2,  and  e*  have  a  multi¬ 
variate  normal  distribution  with  mean 
vector  0  and  covariance  matrix 


Define  the  dummy  variable 
Di  -  1  if  IL  >  0 

-  0  if  <  0 


(1) 

(2) 


-1 


Without  loss  of  generality,  let  us  con¬ 
sider  the  first  equation  in  the  two  re¬ 
gimes  and  assume  thac  the  endogenous  var¬ 
iables  occurring  in  both  are  the  same. 

Let  the  observed  endogenous  variables  be 


equation  (10)  can  now  be  written  as 

yli  “  °iy2i  {8112  -  S212,+  •" 

+  Diypi  (eilp  '  S21P)  +  y21S212  +  ••• 


Let  us  write  the 
two  regimes  as 


if  -  1 
if  Di  =  0 

first  equation  in  the 


(7) 


ylli  *  8112y12i  +  S113y13i  + 
+  6llpylpi  +  YllZi  +  Elli 

and 


y21i  =  8212y22i  +  8213y23i  +  "* 

+  ®21py2pi  +  Y21zi  +  e21i  <9 

These  two  equations  can  be  combined  into 
a  single  equation  and  written  as: 


+  ypi821p  +  Vi<yll'y2l)  +  V21 

+  ^i  ^°21*  ~  °11*^  +  ^li 
where  E(§^)  -  0. 

3.  Estimation  Under  D-Methods 


To  estimate  equation  (12)  and  develop 
valid  statistical  tests  we  will  extend 
the  two-stage  estimations  method  that 
Heckman  (197§a)  suggested.  First,  we  get 
an  estimate  6  of  4  by  probit  maximum  like¬ 
lihood  applied  to  equation  (4) . 

Second,  we  estimate  the  reduced  form  equa¬ 
tions  for  the  endogenous  variables  in  (8) 
and  (9)  by  the  following  two  stage 
method. 


yli  ■  Diy2i(8112  -  ®212>  +  •" 

The  reduced  form  equation  for  y 
can  be  written  as 

+  Diy£i  (8llp  -  W 

yiji  *  Vij  +  viji 

+  y2i8212  +  +  ypi821p 

4, 

Since  E(Vxji  |  Ij  >  0)  ♦  ■  °i;j.  -*7 

+  Diz.(yn  -  y21)  +  z.y21 

where  =  Cov(v1;j,E*) 

Iji 

(13) 


*  Dielli  +  <1-Di)  e21i 


(10) 


We  can  write  (13)  as 


Let  us  define 


♦i  »  «(x|5) 

and 

m  »(x!J| 


(ID 


where  xi  is  the  ith  observation  of 

x  0-  the  standard  normal  density 

i=l,  2..N  1 

evaluated  at  x|6  and  is  the  cumulative 
normal : 


♦i 


_ 1_ 

/TiT 


du 


Note  that  |  >  0)-  - 

^i 

where  °ell*  *  Cov*en'  e*) 

Similarly  E(e2U  I  lL  1  0)  -  0e21* 

where  o21*  ■  Cov(e21,  e*) 

Hence,  the  composite  residual  in  (10)  has 
the  expectation 

E,DiClU  *  (1-Di>  SliJ  “ 

<ac21»  "  °ell*)  *i 


Diylji  *  ♦iZiwlj  ~  °lj*  *i  +  nlji 

(14) 

where  EJnj^)  *  0. 

We  estimate  equation  (J4)  by  qrdinary 
least  squares  substituting  4^  and  4^  for 
4^  and  4^  respectively. 

Here  4^  and  4^  are  the  same  as  4^ 
and  4^  with  4  substituted  for  6.  This 

two-stage  estimation  gives  us  consistent 
estimates  and  and  we  can  con¬ 

struct  the  instrumental  variable 

Diylji  ”  *izi*lj  "  °ij**i  ll5) 

From  the  reduced  form  equations  for  y, 
and  y2jA,  one  has 

yji  -  ♦iVlj  +  <l“*i»Vaj  + 

(02j*-°lj»^i  +  vji 

where  E (Vj ^ )  »  0  (16) 

We  estimate  equation  (16)  by  ordinary 
least  squares  substituting  4^  and  $,  for 
4j,  and  4^  to  get  consistent  estimates  S2j, 

i?2j  and  02j«  “  o^j*. 
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*T 

f 


How  we  construct  instrumental  vari¬ 


ables  for  yV  in  (12)  as  follows 


y5i»Vi’Tij+(l--i)*i’'2j+'i»,2j-lj*)(17) 

With  these  instrumental  variables  con¬ 
structed,  one  can  estimate  the  structural 
coefficients  from  the  following  modified 
equations : 


»ii 


Diy2i<6112_8212)  + 


+  Diypi(6llp-621p’ 


+ 


2i“212  +  •••  +  ypiS21p 


+  ^i2i(YirY21)  +  ZiY21 


(18) 


+  *i(0e2l*‘°^i*)  +  ni 


'ell*'  ’  ''li 
where  nu  -  5n  ♦  (D^*.)  *i(Yu-Y21> 


+  (V!i»  zi<Yn“Y2i> 


+  <W(o*2r°ell*) 


(19) 


and 


’li 


DiElli  +  (1-Di)e21i 


"  (ae21*“aell*,<,>i 
In  matrix  notation,  let 
W  =  [Dy*, Dy *,..., Dy£, 
yj, . . .y*,*z,  z , S] 


(20) 


and  8  (B^^-  ^212 f  ‘  ^ 


(21) 

Ilp-P21p*  *'212  * 


'821p'Yll”Y21' 


Y21'ae21*  "  °ell*) 


»  A  A  A  ‘ 

Also  let  W  =  [Dy*,Dy$, . . . ,Dy*,y$ . 


y*,«Z,  Z,  ♦) 


Then,  the  instrumental  variable  estima¬ 
tor  of  8  is 


It  can  be  shown  that  in  this  case 
given  the  way  the  instrumental  variables 
are  constructed  as  defined  by  equations 
(17)  and  (15) ,  the  instrumental  variable 
estimator  is  the  same  as  the  two-stage 
squares  estimator  so  that 

6  «  (W’W)"1  W'y*  . 

^  (*za) 


The  instrumental  variable  estimator 
can  be  easily  shown  to  be  consistent. 

The  asymptotic  covariance  matrix  of  this 
estimator  is  derived  in  the  appendix.  It 
is  the  last  expression  in  the  appendix. 
Since  it  involves  the  definitions  of  sev¬ 
eral  expressions,  it  is  not  repeated  here. 


4 .  An  Empirical  Illustration 


The  example  we  consider  is  the  esti¬ 
mation  of  the  effects  of  college  educa¬ 
tion  on  earnings  based  on  the  project 


TALENT  data.  More  detailed  discussion  of 
the  data  and  their  use  in  this  problem 
can  be  found  in  Kenny  (1977)  and  Kenny  et. 
al.  (1979).  We  will  here  summarize  all 
the  salient  features  of  the  data  and  the 
model  relevant  for  our  discussion  of  the 
D- method . 

Let  S  =  years  of  college  education 
E  =  earnings 

X,  Z-y,  Zj  are  sets  of  (possibly 

overlapping)  exogenous  variables 
defined  in  Table  1. 

The  switching  simultaneous  equation  sys¬ 
tem  we  consider  is  the  following: 

I  =  X<5  -  e* 

5  =  Zl6l  = 

iff  I  >  0 

E  =  YS  +  Z2B2  +  e12 

S  *  0 

P  _  -  D  .  .  otherwise 

E  Z2®3  +  e22 

I  is  an  indicator  function  such  that  if 
I  >  0,  the  individual  goes  to  college 
and  if  I  <  0  the  individual  does  not  go 
to  college. 

The  variables  in  the  choice  function 
and  the  schooling  equation,  i.e.,  vari¬ 
ables  in  the  sets  X  and  Z.  are:  Rural, 
Split,  MATH,  MOCWAG,  EDMAE,  EDFEM,  CHIL. 
(defined  in  Table  1) . 

The  variables  affecting  the  earnings 
functions,  i.e.,  the  variables  in  the 
set  Z2  are:  Rural,  MATH,  Marital  Status 
and  HC  (defined  in  Table  1) . 

The  variable  S  has  been  coded  as  fol¬ 
lows: 

0  for  those  with  no  college  education 

4.0  for  Bachelor's  degree 

5.5  for  Master’s  degree 

6.0  for  six  years  college 

8.0  for  Ph.D. 

with  numbers  between  0  and  4  for  those 
without  bachelor's  degree  depending  on 
when  they  dropped  out  of  the  4  year  pro¬ 
gram. 

We  will  briefly  discuss  the  reason¬ 
ing  behind  the  choice  of  the  variables 
in  the  sets  X,  Z^,  Z2< 

First,  regarding  the  variables  X  and 
Zj  affecting  the  choice  of  whether  or  not 
to  go  to  college  and  the  years  spent  in 
college,  high  school  students  coming  from 
rural  areas  are  less  likely  to  go  to  col¬ 
lege  or  spend  more  years  in  college. 

Hence  we  included  the  'Rural'  as  a  dummy 
variable.  If  children  are  viewed  as  a 


(23) 

(24) 

(25) 


(26) 
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good  whose  returns  are  specific  to  mar¬ 
riage,  then  divorce  would  lower  the  de¬ 
mand  for  child's  educational  attainment. 
Hence  we  included  the  variable  'Split* 
as  an  explanatory  variable  and  its  coeffi¬ 
cient  is  expected  to  be  negative.  The 
variable  MATH  measures  both  inherited  and 
achieved  cognitive  skills,  hose  with 
higher  MATH  scores  may  learn  more  rapidly 
than  others  and  hence  m^y  find  it  optimal 
to  get  a  higher  degree.  The  variables 
MOCWAG,  EDMAL  and  EOFEM  need  no  justifi¬ 
cation.  We  included  the  variable  CHIL 
because  an  increase  in  the  number  of 
children  in  the  parents'  household 
increases  the  total  cost  of  schooling  for 
the  family  leading  to  a  reduction  in  the 
number  of  years  spent  in  school. 

As  for  the  variables  Z^  affecting 

the  earnings  functions,  the  variable  Rural 
is  a  crude  proxy  for  cost  of  living. 
Individuals  who  grew  up  in  rural  areas, 
where  cost  of  living  (and  thus  wages)  is 
low,  are  more  likely  to  live  in  rural 
areas  than  others.  The  variable  MATH 
measures  cognitive  achievement  and  hence 
should  be  positively  correlated  with 
wages.  The  variable  Marital  Status  is 
expected  to  be  positively  correlated 
with  wages  if  we  argue  that  married 
people  have  greater  responsibilities  and 
hence  are  likely  to  work  harder.  Another 
argument  is  that  the  marriage  market  may 
recognize  that  some  men  are  better  'pros¬ 
pects'  (have  higher  expected  wages)  than 
others.  Finally  the  variable  HC  is  used 
as  a  proxy  for  on  the  job  experience. 

As  the  number  of  jobs  held  (i.e.,  HC) 
increases,  the  amount  spent  with  the  cur¬ 
rent  firm  is  likely  to  decrease.  Hence 
HC  is  expected  to  be  negatively  corre¬ 
lated  with  earnings. 

Table  2  presents  estimates  of  the 
parameters  &  in  equation  (23) .  All  the 
coefficients  have  the  expected  signs. 

Table  3  presents  the  maximum  likelihood 
estimates  of  the  parameters  in  equations 
(24),  (25)  and  (26).  Table  4  presents 
estimates  of  the  parameters  of  the  earn¬ 
ings  functions  obtained  by  the  D-method 
described  in  this  paper.  The  standard 
errors  have  been  computed  from  the  formula 
in  the  appendix. 

Our  objective  in  using  the  D-method 
is  to  test  which  coefficients  in  the 
earnings  functions  are  significantly  dif¬ 
ferent  between  the  two  groups  -  those 
with  college  education  and  those  without. 
The  equations  in  Table  4  suggest  that  it 
is  only  variables  HC  and  MATH  that  are 
significantly  different;  'Marital  Status' 
and  'Rural*  are  not. 

The  example  we  have  used  here  as  an 
illustration  for  the  D-method  is  one 
where  maximum  likelihood  estimation  is 
feasible.  There  would,  however,  be  many 


practical  instances  where  the  use  of 
maximum  likelihood  procedures  would  be 
ruled  out  because  of  the  large  size  of 
the  model  (one  example  is  the  model  con¬ 
sidered  by  Nelson  and  Olson  (1978]).  In 
such  cases  the  D-method  suggested  in  this 
paper  would  be  useful  to  apply  tests  of 
significance  for  equality  of  coefficients 
between  the  two  regimes. 

5.  Conclusions 

t: 

The  paper  presents  an  extension  of 
the  Goldfeld-Quant  D-method  to  switching 
simultaneous  equation  systems  with  endo¬ 
genous  switching.  The  method  is  useful 
to  apply  tests  for  equality  of  coeffi¬ 
cients  across  regimes  in  switching  simul¬ 
taneous  equation  systems  with  selectivity 
The  paper  derives  the  valid  instrumental 
variables  to  use  in  this  situation,  the 
asymptotic  covariance  matrix  of  the  sug¬ 
gested  estimator  and  gives  an  empirical 
illustration. 
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APPENDIX 

In  this  appendix  we  will  derive  the 
covariance  matrix  of  the  two-stage  least 
squares  estimator  (22a) . 


We  have 

(0-0) 


(W’W)-1  w'n 


where  n^,  the  ith  component  of  is 

given  by  equation  (19)  .  Let  us  first  de¬ 
fine  the  N  X  N  diagonal  matrices 

B  =  Diag  d') 

C  =  Diag  (♦i , , *  *  - (2') 
F  =  Diag  [x^6,x^6, . . ,x^6]  ( 3 * ) 

H  =  Diag  [Zjfr ll-^21 ^ ' ' * ‘ ^?Al 1~^21 ) ^ (^ ' ) 


V  =  Diag  [4.<l-4.  ),...4H(1-4M)1 


♦i 


*N 


Diag  [*1(1-^1)",'*n(1-4n)1 
1=  Diag  1^(1-^)  '•••  ^U^)] 
Actually  C2V-1  and/s. 


and 

= 


(S') 

(6') 

(7') 


v  .  miu  '  ' .  —  CV  bUt 

we  have  defined  these  two  diagonal  matrices 
separately  for  convenience. 

Also  let  X  be  the  N  X  L  matrix  of 
observations  on  x, ,  i.e. 


X  = 


N 

Finally  define 
G  -  H-[H+(oell* 

CX(XV'0C)-1X'/'S 


•°e21*)Fl 


(8') 


(9') 


Now  from  equation  (19),  we  have 


o  =  (w'w)-1  w'nj^ 

=  (W'W)-1  W'  l5j+H(D-4)-H(4->) 

+  (°ell*  -  °e2l*)^-'>)l  <10,> 

where  5^,  D,  4,  4,  4,  are  all  column 
vectors  with  their  i-th  elements  5^,  D^, 
*i'  *i'  9i»  ♦i  respectively. 

Expanding  4^  and  4^  up  to  the  first 
power  of  (6-6)  we  get 

34 

♦i  *  +[yyil'  («-«)=4i+4ix|  (6-6) 

Hence  (4-4)=  CX  (6-6)  (11*) 

Similarly 

4  - 4  =  -  FCX  (6-6)  (12') 

Hence  we  have  2  from  equation  (10‘) 

0  -0  -  (W'W)  L  (5.  +  H (D-4) 

*  d3') 

-  HCX  (6-5)-(0Ell*-ae21<[)FCX(6-6)  J 

Now  what  we  will  do  is  express  (6-6)  in 
terms  of  (D-4)  so  that  all  we  have  to  get 
is  Var  (Sj),  Var  (D-4)  and  Cov(5^,D-4). 

Following  Ameniya  (1977)  we  can  write 

(6-6)  Nx'/\X)-1  X’/^1(D-4)  (14') 

Substituting  this  in  (13 ’)  we  have 

(0-0)*(W'W)-1W'  [51+G(D-4)  ]  (15') 

where  G  is  defined  in  ( 9 ' ) . 


From  equation  (20)  in  the  text,  we 

have 

5li  =  e21i+Di  (elli-e21i) 

"  (ae21*-0ell*) ^i  • 

Hence  Var  ( 5 1:i )  =  E  ( § ) 

=  a221+(a211-o221)4. 

+  ‘°e221*-0ell.)(xi6)*i 
°c21*-0ell*)  *i 


Or  Var  (5^  -  a2^  +  (o^-o^)  B 

+  (oc21*“  °ell*)  FC 
-  <0C21*  -  <16’> 


Var  (D-4)  -  v  where  V  is  defined  in  (S') 
and  hence  Var  [G(D-4)J  -  GVG'  (17') 
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Finally,  E(51(D-«)'J=  -°£u*  c 

+  <°cU.  "  °e2l*^  80  <18'> 

Hence,  collecting  the  expressions  (16*), 
(17 • ) ,  and  (18')  we  have 


Var  (6-0)  =  (W'W)-1W' (o^21  IN 

+  (0?21  XN  +  '"liltin'  8 
+  <0?21*-o211*>  FC 


‘0e21*-°£ll*)2  c2  +  GVG' 


-  0£ll*  (CG’+GC) 

+  (°E11*-0£21*) (BCG'+GCB)] 

W  (W'W)-1 

where  the  expressions  B,  C,  F,  G  etc.  are 
all  defined  in  equation  (l1)  to  (8'). 

FOOTNOTES: 

1Richard  Quandt  kindly  pointed  out  to  us 
that  there  is  an  algebraic  error  in  equa¬ 
tion  (3.8),  p.  425  of  the  Maddala-Nelson 
paper.  The  two  denominators  in  4  should 

bG  (y22°22  +  °32>1/2  and  (*!\2+032,1/2 
instead  of  (°22+a32)1/2  and  (°12+°32)l/2 
respectively. 

2 

In  the  discussion  that  follows,  the  nota¬ 
tion  A  means  that  the  two  expressions  have 

the  asymptotic  distributions. 

TABLE  1 

LIST  OF  VARIABLES 

E  =  Natural  logarithm  of  the  hourly 

wage  (in  cents)  in  1971 
S  =  Years  of  college  education  in  1971 

MATH  =  Score  on  a  composite  of  mathe¬ 
matics  achievement  tests  in  1960 
SPLIT  *  0  if  children  living  with  both 
mother  and  father  in  1960 
■  1  Otherwise 

EDMAL  »  Years  of  father's  education 
EDFEM  ■  Years  of  mother's  education 
CHIL  *  Number  of  living  children  (in 

1960)  in  the  family  in  which  the 
respondent  grew  up 

RURAL  *  1  if  pupils  in  grades  9-12  came 
from  an  area  primarily  small 
town  (under  5,000  people)  or 
rural  farm 
»  0  Otherwise. 


TABLE  1  (cont.) 

MOCWAG  “  (mean  occupational  wage)  Log  of 
mean  earnings  (as  of  1960)  of 
full-time  workers  in  father's 
occupation 

HC  =  Number  of  jobs  held  between  1965 
and  1970 

N  =  The  total  sample  size  ■  1373 

N,  =  Number  of  observations  in  the 

1  first  group  (those  who  went  to 

college  =  909) 

N,  =  Number  of  observations  in  the 
*  second  group  (those  that  did  not 

go  to  college  ■  464) 

TABLE  2 

PROBIT  EQUATION  TO  EXPLAIN  THE 
DECISION  OF  WHETHER  OR  NOT  TO  GO 
TO  COLLEGE 


Variable 

Coeff.  Est. 

S.E. 

Constant 

-5 

.259 

2.108 

Rural  (Dummy) 

.164 

.093 

Split 

- 

.078 

.145 

MATH 

.015 

.001 

MOCWAG 

.922 

.589 

EDMAL 

.044 

.017 

EDFEM 

.029 

.019 

CHIL 

.057 

.021 

TABLE  4 

EARNINGS  FUNCTIONS  WITH  D-METHOD 

Variables 

Equation 

1  Equation  2  Equation 

Constant 

5.72467 

5.69896 

5.70961 

(0.12918) 

(0.08953) 

(0.09773) 

Marital 

0.14620 

0.15656 

0.15511 

Status 

(0.06281) 

(0.01298) 

(0.01808) 

Rural 

-0.07161 

-0.05352 

-0.05030 

(0.07239) 

(0.02185) 

(0.02416) 

MATH 

-0.00359 

-0.00377 

-0.00336 

(0.00231) 

(0.00206) 

(0.00193) 

HC 

0.03196 

0.03184 

0.03234 

(0.01304) 

(0.01275) 

(0.01315) 

♦ 

0.74827 

0.82369 

0.66080 

(0.48986) 

(0.35637) 

(0.20604) 

4 

-0.12980 

-0.15768 

(0.34249) 

(0.27731) 

4*Marital 

0.01454 

(0.08599) 

4  *Rural 

0.03258 

(0.10696) 

4  MATH 

0.00560 

0.00589 

0.00507 

(0.00306) 

(0.00262) 

(0.00226) 

4*HC 

-0.06444 

-0.06335 

-0.06671 

(0.02026) 

(0.01920) 

(0.01894) 

S 

0.05141 

0.05774 

0.04202 

(0.04710) 

(0.03536) 

(0.02590) 

(Table  3 

is  available  from  the 

authors 

upon  request.) 
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